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B —> 7 t££ form factors for new-physics searches from lattice QCD 


Jon A. Bailey, 1 A. Bazavov, 2 C. Bernard, 3 C.M. Bouchard, 4 C. DeTar, 5 Daping Du, 6, * 

A.X. El-Khadra, 7 E.D. Freeland, 8 E. Gamiz, 9 Steven Gottlieb, 10 U.M. Heller, 11 A.S. Kronfeld, 12,13 
J. Laiho, 6 L. Levkova, 5 Yuzhi Liu, 14 E. Lunglii, 10, ' P.B. Mackenzie, 12 Y. Meurice, 2 E. Neil, 14,15 
Si-Wei Qiu, 5 J.N. Simone, 12 R. Sugar, 16 D. Toussaint, 17 R.S. Van de Water, 12 A and Ran Zhou 12 

(Fermilab Lattice and MILC Collaborations) 

1 Department of Physics and Astronomy, Seoul National University, Seoul, South Korea 
2 Department of Physics and Astronomy, University of Iowa, Iowa City, IA, USA 
3 Department of Physics, Washington University, St. Louis, MO, USA 
4 Department of Physics, The Ohio State University, Columbus, OH, USA 
5 Department of Physics and Astronomy, University of Utah, Salt Lake City, UT, USA 
6 Department of Physics, Syracuse University, Syracuse, NY, USA 
' Department of Physics, University of Illinois, Urbana, IL, USA 
8 Liberal Arts Department, School of the Art Institute of Chicago, Chicago, IL, USA 
9 CAFPE and Departamento de Fisica Tedrica y del Cosmos, Universidad de Granada, Granada, Spain 
10 Department of Physics, Indiana University, Bloomington, IN, USA 
11 American Physical Society, Ridge, NY, USA 
12 Fermi National Accelerator Laboratory, Batavia, IL, USA 
13 Institute for Advanced Study, Technische Universitat Miinchen, Garching, Germany 
14 Department of Physics, University of Colorado, Boulder, CO, USA 
15 RIKEN-BNL Research Center, Brookhaven National Laboratory, Upton, NY, USA 
16 Department of Physics, University of California, Santa Barbara, CA, USA 
11 Physics Department, University of Arizona, Tucson, AZ, USA 
(Dated: August 14, 2015) 

The rare decay B —> n£ + £~ arises from b —> d flavor-changing neutral currents and could be 
sensitive to physics beyond the Standard Model. Here, we present the first ab-initio QCD calculation 
of the B —> 7r tensor form factor fr ■ Together with the vector and scalar form factors /+ and fo 
from our companion work [J. A. Bailey et al, Phys. Rev. D 92 , 014024 (2015)], these parameterize 
the hadronic contribution to B —> n semileptonic decays in any extension of the Standard Model. 

We obtain the total branching ratio BR(R + —> n + p + p.~) = 20.4(2.1) x 10 -9 in the Standard Model, 
which is the most precise theoretical determination to date, and agrees with the recent measurement 
from the LHCb experiment [R. Aaij et al., JHEP 1212 , 125 (2012)]. Note added: after this paper 
was submitted for publication, LHCb announced a new measurement of the differential decay rate for 
this process [T. Tekampe, talk at DPF 2015], which we now compare to the shape and normalization 
of the Standard-Model prediction. 

PACS numbers: 13.20.He, 12.38.Gc, 12.15.Mm 


Motivation — Hadron decays that proceed through 
flavor-changing neutral currents may be sensitive to new 
physics, because their leading Standard-Model contribu¬ 
tions are loop suppressed. Here we study the semilep¬ 
tonic decay B —» 7 t£ + 1~, which proceeds through a 
b -A d transition. Hadronic effects in this decay are 
parametrized by three form factors. In this Letter, we 
present the first ab-initio QCD calculation of the ten¬ 
sor form factor fo, based on lattice-QCD work that 
also yielded the vector and scalar form factors, f + and 
fo [1]. Lattice QCD has several advantages over other 
approaches to the form factors [2-9], particularly in pro¬ 
viding a path to controlled uncertainties that can be sys¬ 
tematically reduced [10]. 

The LHCb experiment recently made the first obser¬ 
vation of B + —> 7 t + pL + [11], while the B -factories 
have set limits on the e + e~ and t + t~ channels [12-14], 
Below we present the first calculations of B —> ir£ + £~ 
{£ = e,/x,r) observables in the Standard Model using 


form factors with fully controlled uncertainties. 

The form factors /+, fo, and fr suffice to parame¬ 
terize B —> 7r decays in all extensions of the Standard 
Model. New physics from heavy particles—such as those 
appearing in models with supersymmetry [3, 15-17], a 
fourth generation [18], or extended [15, 19-23] or com¬ 
posite [24] Higgs sectors—alter Wilson coefficients in the 
effective Hamiltonian pertaining to particle physics below 
the electroweak scale [25-28]. Whatever these unknown 
particles may be, the hadronic physics remains the same. 

Lattice-QCD calculation — Our work on fr(<l 2 ) was 
carried out in parallel with f + (q 2 ) and /o(? 2 )- Our aim 
in Ref. [1] was a precise determination of the Cabibbo- 
Kobayashi-Maskawa (CKM) element V u i ,|, and every 
step of the analysis was subjected to many tests. Fur¬ 
ther, two of the authors applied a multiplicative offset 
to the numerical data at an early stage. This “blinding” 
factor was disclosed to the others only after finalizing the 
error analysis. Full details of the simulation parameters, 
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analysis, and cross-checks are given in Ref. [1]. 

Our calculation uses ensembles of lattice gauge-field 
configurations [29] from the MILC Collaboration [30-32] , 
which are generated with a realistic sea of up, down, and 
strange quarks. In practice, the up and down sea quarks 
have the same mass, and the strange-quark mass is tuned 
close to its physical value. The statistics are high, with 
600-2200 gauge-field configurations per ensemble. The 
physical volume is large enough that we can repeat the 
calculation in different parts of the lattice, thereby qua¬ 
drupling the statistics. We use four lattice spacings rang¬ 
ing from 0.12 fm to 0.045 fm to control the extrapolation 
to zero lattice spacing. 

The tensor form factor is defined via the matrix ele¬ 
ment of the b —> d tensor current ida^b: 

(n{p w )\ida^b\B(p B )) = Mq 2 ), (1) 

where p B and p w are the particles’ momenta and q = 
Pb —p-w is the momentum carried off by the leptons. The 
Lorentz invariant q 2 is related to the pion energy in the 
B-meson rest frame via E n = (M| + M 2 — q 2 )/2M B . In 
the finite volume that can be simulated on a computer, 
E n takes discrete values, dictated by the spatial mo¬ 
menta p v compatible with periodic boundary conditions. 
Because statistical and discretization errors increase with 
pion momentum, we restrict \p^\ < |27r(l, 1,1)/L|. The 
resulting simulation range of E n < 1 GeV is signifi¬ 
cantly smaller than the kinematically allowed range of 
E n < 2.5 GeV. Extending this discrete set of calcula¬ 
tions into the full q 2 dependence is the central challenge 
of this work, and is met in two steps. 

The two light quarks (up and down) have a mass larger 
than it should be, but the range simulated is wide and 
the smallest pion mass is 175 MeV, close to Nature’s 
140 MeV. Therefore, we can apply an effective field the¬ 
ory of pions—chiral perturbation theory—to extrapolate 
the simulation data to the physical point. We use a form 
of chiral perturbation theory adapted to lattice QCD, 
with additional terms describing the lattice-spacing de¬ 
pendence [33, 34] and with modifications needed for en¬ 
ergetic final-state pions [35]. As discussed in Ref. [1], we 
try several fit variations. For example, we replace the 
loop integrals with momentum sums appropriate to the 
finite volume, finding negligible changes in the results. 
Our final fit includes next-to-next-to-next-to-leading or¬ 
der analytic terms and terms to model the discretization 
errors of the heavy quark. The latter come from an ef¬ 
fective field theory for heavy b quarks [36-38]. 

Figure 1 shows the q 2 dependence of the errors after the 
chiral-continuum extrapolation just described. Table I 
gives a numerical error budget for /t(? 2 = 20 GeV 2 ). 
The largest uncertainty comes from the statistical errors, 
as increased during the chiral-continuum extrapolation. 
This error is under good control for q 2 corresponding to 
the spatial momenta that we simulate, but grows large 


TABLE I. Error budget in per cent for fx(q 2 = 20 GeV 2 ). 
The first error incorporates statistical errors from the simu¬ 
lation and systematics associated with the chiral-continuum 
fit. The last column emphasizes how the error varies with q 2 . 


Source of error 

Sfx 

q 2 dependence 

Statistics © yPT © HQ ffi g w 

3.8 

important 

Scale ri 

0.5 

negligible 

Nonperturbative matching Z v a , Z v 4 

0.7 

negligible 

Perturbative matching px 

2.0 

none 

Heavy-quark mass tuning Kb 

0.4 

none 

Light-quark mass tuning mi , m s 

0.5 

negligible 

Total ( Quadrature sum of above) 

4.4 

important 


elsewhere. 

The subdominant errors are as follows. To convert 
from lattice units to physical units, we introduce a 
physical distance ri, which is defined via the force be¬ 
tween static quarks [39, 40]. We use it to form physi¬ 
cal, dimensionless quantities, which are the input data 
for the chiral-continuum fit. At the end, we set r\ = 
0.3117 ± 0.0022 fm [41] based on a related lattice-QCD 
calculation of rif n [42] and the pion decay constant 
f n = 130.41 MeV [43]. To propagate the parametric 
uncertainty in r\ to fr, we repeat the fit shifting ri by 
±1 a ri , leading to the second line in Table I. 

In lattice gauge theory, the tensor current does not 
have the normalization used in QCD phenomenology. We 
obtain most of the normalization nonperturbatively [44] 
from b —> b and d —> d transitions with the vector current, 
with statistical errors below 1%. Another matching fac¬ 
tor pt remains, but, by design and in practice, it is close 
to unity. We calculate p B at the renormalization scale 
p = m^ p 0 ie through first order in the QCD coupling a s . 



17 18 19 20 21 22 23 24 25 26 

q 2 (GeV 2 ) 


FIG. 1. (color online) Error budget for fx as a function of q 2 
for the range of simulated lattice momenta. The filled bands 
show the relative size of each error contribution to the total. 
The quadrature sum is shown on the left y axis and the error 
itself, in per cent, on the right. 
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We estimate the resulting error of order a 2 after remov¬ 
ing a logarithmic dependence on the matching scale p, 
which is present in continuum QCD too. We then exam¬ 
ine how the one-loop coefficient depends on heavy-quark 
mass, identifying the largest value, Pj^ max - Finally, we 

estimate the error in px to be 2a a |p^ max |, evaluating a s 
on the second-finest lattice with a « 0.06 fm. This yields 
the 2% perturbative-matching uncertainty in Table I. 

The last two uncertainties arise as follows. When gen¬ 
erating data, we choose the simulation quark masses 
based on short runs and previous experience. The full 
analysis yields better estimates. To correct the simu¬ 
lation 6-quark mass a posteriori , we recompute fx on 
one ensemble with two additional values of the bare 6- 
quark mass. Using the slope from all three mass values, 
we interpolate the data for fx slightly from the produc¬ 
tion 6-quark mass to the physical value. This leaves an 
error due to the uncertainty in the size of the 6-quark 
mass correction. The details for fx are nearly identical 
to those for f + [1], leading to the same estimate, 0.4%, 
for this error. The light-quark mass dependence is em¬ 
bedded in the chiral-continuum extrapolation, described 
above. The parametric uncertainty from the input light- 
quark mass [32] is propagated to fx by repeating the fit 
with ±lcr m shifts to these parameters, and is given in 
the penultimate line of Table I. 

The final line in Table I and the upper edge of the stack 
in Fig. 1 represent the quadrature sum of the systematic 
uncertainties with the chiral-continuum fit error. 

Extension to all q 2 — To extend fx in the chiral- 
continuum limit from the range of simulated lattice mo¬ 
menta to the full kinematic range, 0 < q 2 < {Mb — M n ) 2 , 
with controlled errors, we use a method based on the an¬ 
alytic structure of the form factor. 

In the complex q 2 plane, fx{q 2 ) has a cut for timelike 
q 2 > t + = {Mb + M n ) 2 and a pole at q 2 = Mg* but is 
analytic elsewhere. The variable 


z{q 2 ,t 0 ) 


Vt + - q 2 - y/t+ -t 0 

\A+ ~ q 2 + y/t + -t 0 


( 2 ) 


maps the whole q 2 plane into the unit disk, with the 
cut mapped to the boundary and the semileptonic region 
mapped to an interval on the real axis. Unitarity bounds 
then guarantee that an expansion of fx in 2 (with the B* 
pole removed) converges for |z| < 1 [45-48]. Following 
Bourrely, Caprini, and Lellouch (BCL) [49], we factor out 
the B* pole and expand in 2 : 


(1 -q 2 /M 2 *)f T {z) 


N z -1 

n—0 


'W* 


N z 


______ ( 3 ) 

choosing t 0 = {M B +M V ){y/M B — \/ M^) 2 to minimize 1 2 | 
in the semileptonic region. Although Eq. (3) was derived 
for the vector form factor f + , we use it for the tensor form 


TABLE II. Best-fit values b n with total errors and correlation 
matrix p nm of the N z = 4 BCL 2 expansion of fr ■ The lower 
two panels show correlations with the /+ and /o coefficients 
in Table XIV of Ref. [1] obtained from ab-initio QCD. 


Fit: 

0.393(17) 

-0.65(23) 

—0.6(1.5) 0.1(2.8) 

P 

bo 

bi 

bi 

63 

bo 

1.000 

0.400 

0.204 

0.166 

b'f 


1.000 

0.862 

0.806 

hi 



1.000 

0.989 

63 




1.000 

&o 

0.638 

0.321 

0.123 

0.084 

bi 

0.321 

0.397 

0.162 

0.109 

bi 

0.114 

0.202 

0.198 

0.179 

bi 

0.070 

0.152 

0.192 

0.180 

&o 

0.331 

0.136 

0.089 

0.073 

6? 

0.203 

0.313 

0.198 

0.162 

b°2 

0.204 

0.268 

0.186 

0.155 

63 

0.151 

0.203 

0.169 

0.149 


factor fx because the two form factors are proportional 
to each other at leading order in the 1/mb expansion. 

We determine the with a functional method 
connecting the independent functions of the chiral- 
continuum fit with the first several powers of 2 [1]. Our 
preferred fit uses N z = 4; adding higher-order terms in 
2 does not significantly change the central value. Ta¬ 
ble II presents our final result for fx as coefficients of 
the N z = 4 BCL 2 fit and the correlation matrix be¬ 
tween them, where the errors include statistical and all 
systematic uncertainties. This information can be used 
to reconstruct fx{q 2 ) over the full kinematic range. Ta¬ 
ble II also provides the (mostly statistical) correlations 
between fx, /+, and / 0 . Figure 2 shows the extrapola¬ 
tion of fx to q 2 = 0. Table II and Fig. 2 represent the 
first main result of this Letter. 

Implications — The largest contribution in the Stan¬ 
dard Model to the amplitude for B —► 7 tI + I~ is pro¬ 
portional to the vector form factor. Assuming that new 
physics does not contribute significantly to the tree-level 



FIG. 2. Ab-initio result for fr{q 2 ) from lattice QCD. 



















4 


TABLE III. Correlations between BCL coefficients for fx with 
those for /+ and /o from Table XIX of Ref. [1] , which include 
experimental shape information from B —> n£v decay. 


p 

60 

b'[ 

62 

63 

bt 

0.514 

0.140 

0.078 

0.065 

bf 

0.111 

0.221 

- 0.010 

-0.049 

bt 

-0.271 

-0.232 

- 0.012 

0.029 

4 

-0.204 

-0.215 

-0.013 

0.023 

b\\ 

0.243 

-0.015 

-0.025 

-0.024 

6 ? 

0.005 

0.134 

0.070 

0.057 

b% 

- 0.002 

-0.034 

-0.032 

-0.030 

b u 3 

-0.044 

-0.061 

0.005 

0.017 


decay B —> 7 t£v, one can use experimental measurements 
of this process to constrain the shape of f + (q 2 ), especially 
at low q 2 . In Ref. [1], we obtain the CKM element \V u b\ 
from a combined z fit to our lattice-QCD results for f + 
and /o and measurements of tbcIT(B —> Tt£v)/dq 2 from 
BaBar [50, 51] and Belle [52, 53]. This joint fit also yields 
the most precise current determinations of /+ and fo■ To 
enable them to be combined with the results for fx from 
Table II, Table III provides the correlations between the 
^-expansion coefficients for all three form factors. The 
correlations are small, because / + contains independent 
experimental information. 

Using fx from this work and /+ and /o just described, 
we show the Standard-Model partial branching fractions 
for B —> 7 t£ + £~ in Fig. 3. Other ingredients are needed 
besides the form factors. We take the Wilson coefficients 
from Ref. [27], the CKM elements from Ref. [55], the me¬ 
son masses and lifetimes from Ref. [43], and the b- and 
c-quark masses from Ref. [7]. To calculate contributions 
that cannot be parameterized by the form factors, we em¬ 
ploy QCD factorization at low q 2 [56-64] and an operator 
product expansion (OPE) in powers of E n /^fq 2 at large 
q 2 [65-72]. Full details will be provided in Ref. [73]. 

Table IV presents numerical predictions for selected 
q 2 bins. The last error in parenthesis contains effects 
of parametric uncertainties in a s , mt , mb, m c \ of miss¬ 
ing power corrections, taking 10% of contributions not 
directly proportional to the form factors; and of vio¬ 
lations of quark-hadron duality, estimated to be 2% at 
high-g 2 [70]. At low g 2 , the uncertainty predominantly 
stems from the form factors; at high q 2 , the CKM ele¬ 
ments |V" t ^Vtb| and form factors each contribute similar 
errors. Figure 3 and Table IV represent the second main 
result of this Letter. 

In the regions q 2 < 1 GeV 2 and 6 GeV 2 < q 2 < 
14 GeV 2 , uu and cc resonances dominate the rate. To 
estimate the total BR, we simply disregard them and in¬ 
terpolate linearly in q 2 between the QCD-factorization 
result at q 2 ~ 8.5 GeV 2 and the OPE result at 
q 2 ss 13 GeV 2 . While this treatment does not yield 
the full branching ratio, it does enable a comparison 
with LHCb’s published result, BR(R + —>■ 7r + /r + /r“) = 



0 5 10 15 20 25 

q 2 [GeV 2 ] 


FIG. 3. (color online) Partial branching fractions for B + —> 
7 r + /r + /r~ (upper panel) and B + —> tv + t + t~ (lower panel) out¬ 
side the resonance regions. Different patterns (colors) show 
the contributions from the main sources of uncertainty; those 
from the remaining sources are too small to be visible. For 
B + —y n + i_i + fj,~, new measurements from LHCb [54], which 
were announced after our paper appeared, are overlaid. 


TABLE IV. Standard-Model predictions for B + —> n + £ + £~ 
partial branching fractions. Those for B° decays can be ob¬ 
tained by multiplying by the lifetime ratio (t b o/t b +)/2 = 


0.463. Errors shown are from the CKM elements, form fac¬ 
tors, variation of the high and low matching scales, and the 
quadrature sum of all other contributions, respectively. 

[^min ? 9max] 

10 a x BR (B+ - 

-+ n + £ + £~) 

(GeV 2 ) 

£ = e, fi 

£ = T 

[0.1, 2.0] 

1.81(11,24,6,2) 


[2.0,4.0] 

1.92(11,22,6,3) 


[4.0,6.0] 

1.91(11,20,6,3) 


[6.0,8.0] 

1.89(11,18,5,3) 


[15,17] 

1.69(10,13,3,5) 

1.11(7,8,2,4) 

[17,19] 

1.52(9,10,2,4) 

1.25(8,8,2,3) 

[19,22] 

1.84(11,11,3,5) 

1.93(12,10,4,5) 

[22, 25] 

1.07(6,6,3,3) 

1.59(10,7,4,4) 


[ 1 , 6 ] 

[15,22] 

4.78(29,54,15,6) 

5.05(30,34,7,15) 4.29(26,25,7,12) 

[4m, 2 , 26.4] 

20.4(1.2,1.6,0.3,0.5) 


23(6) x 10 -9 [11], which was obtained from a similar in¬ 
terpolation over these regions. Our result BR(R + —► 
7 r + /i + /r - ) = 20.4(2.1) x 10” 9 agrees with LHCb, and 
is more precise than the best previous theoretical esti¬ 
mate [7] because we use fx directly, which avoids a large 
uncertainty from varying the matching scale fi. 
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Outlook - The largest uncertainty in our determi¬ 
nation of the B —► 7r form factors is the combined 
error from statistics with chiral-extrapolation and dis¬ 
cretization effects included. We will be able to re¬ 
duce these with calculations on the MILC Collabora¬ 
tion’s recently generated four-flavor ensembles with phys¬ 
ical light-quark masses [74]. LHCb’s measurement of 
BR(R + —> n + ii + fi ~) will improve, and Belle II expects 
to observe the neutral decay mode B° —> tt°£ + £~. If a 
deviation from the Standard Model is observed, our form 
factors can be used to compute other observables such as 
asymmetries, thereby providing information about new 
heavy particles, such as their masses, spin, and couplings. 

Note added: after this paper was submitted for publi¬ 
cation, the LHCb experiment announced a new measure¬ 
ment for the B + —> 7r + /i + /r _ differential decay rate [54], 
The new results are shown in Fig. 3. The large differ¬ 
ence in the lowest q 2 bin is due to the presence of light 
(/?, w, </>) resonances, whose effects are important but can¬ 
not be estimated in a model-independent manner. Given 
the present experimental and theoretical uncertainties, it 
is too early to discern possible new-physics contributions 
to this process. 
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